Strain-induced ferroelectricity in simple rocksalt binary oxides 
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The alkaline earth binary oxides adopt a simple rocksalt structure and form an important family 
of compounds because of their large presence in the earth's mantle [1] and their potential use in mi- 
croelectronic devices [2]. In comparison to the class of multifunctional ferroelectric perovskite oxides 
[3], however, their practical applications remain limited and the emergence of ferroelectricity and 
related functional properties in simple binary oxides seems so unlikely that it was never previously 
considered. Here, we show using first-principles density functional calculations that ferroelectricity 
can be easily induced in simple alkaline earth binary oxides such as barium oxide (BaO) using ap- 
propriate epitaxial strains. Going beyond the fundamental discovery, we show that the functional 
properties (polarization, dielectric constant and piezoelectric response) of such strained binary ox- 
ides are comparable in magnitude to those of typical ferroelectric perovskite oxides, so making them 
of direct interest for applications. Finally, we show that magnetic binary oxides such as EuO, with 
the same rocksalt structure, behave similarly to the alkaline earth oxides, suggesting a route to new 
multiferroics combining ferroelectric and magnetic properties. 



The alkaline earth oxides - MgO, CaO, SrO and BaO 
- have been intensively studied [4] and their electronic 
structures, elasticities, thermal properties and equations 
of state are well established [5, 6]. In view of their rela- 
tively wide band gaps and their compatibility with sili- 
con, they recently aroused some interest as possible gate 
dielectrics to replace Si0 2 in MOSFET devices [2]. Also, 
they have been shown to play a key role as a buffer layer 
in the epitaxial growth of multifunctional perovskite ox- 
ides directly on silicon [7]. As a result, many recent 
works have been devoted to the study of AO oxide thin 
films on silicon substrates, with a special emphasis on 
their lattice mismatch, coherence and band offsets with 
Si [2, 7, 8]. Surprisingly, although many studies of the 
phase diagrams at high pressure have been performed [9], 
the phase diagrams of alkaline earth oxides under epitax- 
ial strain have not been previously reported. 

Ferroelectricity in oxides is usually associated with the 
family of ABO3 perovskites [3]. The ferroelectric phase 
transition in this class of materials is described in terms 
of a displacive transition from a high-symmetry paraelec- 
tric phase by condensation of a polar transverse-optical 
(TO) soft phonon mode below the critical temperature. 
Experimentally, this is characterized by a decrease in fre- 
quency of the TO mode as the ferroelectric transition 
temperature is approached from above; computationally, 
the frequency of the TO mode is found to be imaginary 
in the paraelectric phase, consistent with the absence of 
a restoring force for ionic displacements [10]. As first 
proposed by Cochran [11], the ferroelectric instability 
can be explained from the compensation of short-range 
forces favoring the undistorted paraelectric structure by 
long-range Coulomb interactions favoring the ferroelec- 
tric phase. This was confirmed at the first-principles level 
[12], pointing out that the unusually large destabilizing 
Coulomb interaction yielding the instability in this class 
of compounds is linked to giant anomalous Born effective 



charges, Z* (2.77 e, 7.25 e and -5.71 e for Ba, Ti and O 
in BaTiOa instead of nominal charges of 2.0 e, 4.0 e and - 
2.0 e), themselves produced by the strong sensitivity of O 
2p - metal d hybridizations to atomic displacements [13]. 
The tendency of materials to ferroelectric instability was 
also recently reformulated in the framework of vibronic 
coupling theory [14, 15]. There again, the rearrange- 
ment of electrons through covalent bond formation when 
atoms are displaced, typically resulting in anomalous Z* , 
appears as an important feature to yield ferroelectricity. 

Alkaline earth oxides exhibit anomalous Z*, related to 
O 2p - metal d hybridizations through a mechanism sim- 
ilar to ABO3 compounds [6]: We obtain values of 2.81 e, 
2.49 e and 2.39 e respectively for BaO, SrO and CaO (the 
nominal charge is 2.0 e in all cases). Nevertheless, ten- 
dency to ferroelectricty has never been reported or stud- 
ied in alkaline earth oxides . The ferroelectric instablity 
is known to be strongly sensitive to strain and, recently, 
it was shown that strain engineering can be used to in- 
duce ferroelectricity in otherwise non-ferroelectric ABO3 
compounds [16, 17]. Here we explore whether the same 
strategy can be applied to AO compounds, by studying 
the properties of different binary oxides under epitaxial 
strain using density functional theory calculations (see 
Methods section). 

We begin by calculating the TO phonon frequencies 
of the alkaline earth oxides as a function of epitaxial 
strain to search for the onset of a softening or instabil- 
ity. In their non-distorded cubic rocksalt structure, these 
oxides have a three-fold degenerate TO mode at the T 
point. Under biaxial epitaxial strain, as can be practi- 
cally achieved in thin films on a substrate, the symmetry 
is reduced to tetragonal (space group H/mmm, No. 139), 
inducing a splitting of the cubic TO mode into a single 
A 2 « mode polarized along the [001] tetragonal axis, and 
a two-fold degenerate E u mode polarized either along the 
[100] or [010] axis. 
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We first consider the case of BaO. Figure l.a shows the 
evolution of the A2 U and E u frequencies with respect to 
a wide range of epitaxial strain, r\ = (a s — ao)/ao where 
ao is the relaxed lattice constant of BaO and a s is the 
lattice constant of the substrate. We clearly see that the 
TO mode frequencies are strongly strain-dependent. For 
compressive epitaxial strains (77 < 0), the E M frequency 
increases smoothly whereas the A2 U frequency decreases 
strongly until it becomes imaginary at a critical epitax- 



ial strain i]* ~-1.3% (left vertical line of Figure l.a). A 
similar behavior is observed for tensile epitaxial strains 
(7/ > 0), but here it is the E u frequency that becomes 
imaginary at rff 1 ~1.8% (right vertical line of Figure l.a). 
The existence of a polar TO mode with an imaginary 
frequency is the fingerprint of ferroelectricity: The con- 
densation of either the A2 U or E M unstable mode in BaO 
will induce a spontaneous polarization and yield a ferro- 
electric ground-state. 




FIG. 1: (a) Frequencies of the A21, (red circles) and E u (blue squares) modes of BaO with respect to the epitaxial strain. 
Negative numbers correspond to imaginary frequencies, (b) Energy as a function of polarization for epitaxial BaO. Green 
crosses, blue squares and red circles correspond to an epitaxial strain of respectively -0.6%, -1.63% and -2.67%. (c) BaO 
static dielectric constant eq versus the epitaxial strain (eq x red squares, Eq V green crosses, eg 2 blue circles) and (d) 1/eq versus 
the epitaxial strain (l/so x red squares, I/eq 2 blue circles), (e) Frequencies of the A2 U (red circles) and E u (blue squares) 
modes of BaO/SrO superlattice with respect to the epitaxial strain. Negative numbers correspond to imaginary frequencies, 
(f ) Amplitude of the spontaneous polarization of EuO with respect to the epitaxial strain. Red circles and blue squares are 
respectively for out-of-plane and in-plane orientation of the polarization. The vertical lines of graphs a, c, d, e and f give the 
position of the critical strains rf c and rfc v . 



To determine the ground-state under high epitaxial 
strains, we next performed full atomic relaxations, in- 
tialized by condensing the unstable polar modes, for each 
mistfit strain exceeding the critical value. For compres- 
sive epitaxial strain, the strongest instability corresponds 
to polar distortions along the [001] direction (A2„). The 
condensation of this instability brings the system to a 
new ground-state, hereafter called the c-phase, with a 
non-zero polarization along the [001] direction (P z ). The 
gain of energy with respect to the paraelectric reference 
(AE), the polarization (P) and the tetragonality c/a are 
reported in Table I. At an epitaxial strain of -2.66%, 
these numbers are huge for an otherwise classical dielec- 



tric compound and comparable in amplitude to a typical 
ferroelectric such as BaTi0 3 (P~34 /iCcm -1 , AE~3.8 
meV/atom). Figure l.b shows the evolution of the energy 
with the polarization resulting from the atomic relaxation 
at an epitaxial strain of -2.63% (red circles curve). As im- 
posed by symmetry, the variation of energy is the same 
when the polarization is oriented up or down, so that the 
energy displays the double-well shape characteristic of 
a ferroelectric compound. For comparison, an analogous 
plot is shown at an epitaxial strain of -0.6% (green crosses 
in Figure l.b) where no polar instability is present. Here 
a single-well shape is obtained, indicating that there is 
no energy gain when the compound is polarized. 
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The absence of any other phonon instability in the po- 
lar phase, as explicitly checked from the computation 
of the full phonon dispersion curves, confirms that the 
ground-state indeed corresponds to the ferroelectric c- 
phase when r\ > rf c . As in the case of the ferroelectric 
perovskites, the tetragonality is strongly modified by the 
polarization; we find that the polarization increases the 
tetragonality (c/a=(c/a)o+aP 2 ) of the structure with re- 
spect to the paraelectric reference and the change of c/a 
is continuous at rf c . 
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TABLE I: Amplitude of the spontaneous polarization |P|, gain 
of energy AE with respect to the paraelectric reference, ratio 
between the out-of-plane and the in-plane cell parameter (Ba- 
O distances) of the paraelectric reference (c/a)o and of the 
ground-state c/a for five epitaxial strains where the c-phase 
or the aa-phase are stable. 



Similarly, under tensile epitaxial strains, the conden- 
sation of the unstable E M mode for 77 > rff 1 yields a 
ferroelectric ground-state with an in-plane polarization 
hereafer called the a-phase. The value of the sponta- 
neous polarization, depth of the double well and tetrago- 
nality are reported for the a-phase in Table I. For equiv- 
alent strain amplitude, the spontaneous polarization and 
the depth of the double well are smaller than the cor- 
responding c-phase values, but are still comparable to 
typical perovskite ferroelectrics. The modification of c/a 
with respect to the paraelectric (c/a)o is almost negligi- 
ble with a tendency to decrease when the epitaxial strain 
or P increases (a <0). 

Next, we evaluate the effect of epitaxial strain on the 
functional properties such as the dielectric permittivity 
and piezoelectric response. Figure l.c shows the evolu- 
tion of the three components e xx , £q V and e% z of the static 
dielectric tensor with the epitaxial strain. For ?7=0%, the 
dielectric tensor is isotropic since the structure is cubic. 
For positive and negative epitaxial strain, Eq z becomes 
different from e xx and Eq V . The value of e" (resp. e xx 
and £q V ) diverges around rf c (resp. rj xv ). This divergent 
behavior of £0, exemplified from the linear evolution of 
l/so in Figure l.d, is a usual feature of a displacive fer- 
roelectric phase transition and is related to the softening 
of a transverse optic polar mode (Eq.55 of Ref [18]). 

Since the ground-state is non-centrosymmetric beyond 
the critical epitaxial strain ry c , it will also be piezoelectric. 
In Table II we report the piezoelectric stress constants e^ 



of BaO for selected epitaxial strains in the ferroelectric 
region. For negative epitaxial strains, in the c-phase, the 
dominant coefficient is 633. This coefficient is large and 
is comparable to that of ferroelectric perovskites such as 
PbTi0 3 (3.68 Cm" 2 [19]), BaTiOg (3.64 Cm" 2 [20]) 
or even the widely used Pb(Zr,Ti)0 3 (PZT) alloys (3.4 
C.m~ 2 [21]). In the a-phase for tensile epitaxial strains, 
the largest piezoelectric constants are en, ei2, ei3 and 
e26 • 
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TABLE II: Non-zero piezoelectric stress constants ey 
(Cm -2 ) of BaO at epitaxial strains of -2.67%, -1.63%, 1.99% 
and 3.03%. 



Now that we have demonstrated the possibility of in- 
ducing ferroelectricity in BaO thin films under epitaxial 
strain, we investigate if the phenomenon is specific to 
BaO or a more general feature shared with other alka- 
line earth compounds. Looking at CaO and SrO, we 
find a very similar behavior. Compressive (resp. tensile) 
epitaxial strain makes the A2 U (resp. E u ) mode softer. 
However, since at the bulk level, the cubic TO frequencies 
of SrO and CaO are higher than in BaO, the calculated 
critical epitaxial strains are also larger : rfc ~-5.1% (resp. 
-6.7%) and r%' y ~ 6.1% (resp. 7.7%) in SrO (resp. CaO). 

Although the larger critical strains make CaO and SrO 
less attractive for applications than BaO, the natural 
tendency to ferroelectricity of the class of alkaline earth 
compounds suggests tuning of the ferroelectric properties 
in superlattices made of the repetition of different ox- 
ide layers; such behavior was realized recently in ABO3 
compounds [22, 23]. As a prototypical example, we con- 
sider here a BaO/SrO superlattice built by alternating 
one layer of BaO and one layer of SrO along the [001] 
direction. The lattice mismatch between BaO and SrO 
is relatively large (~7.7%) and will prevent the growth of 
large period superlattices but should not be problematic 
in ultra-short period systems. The value of the relaxed 
in-plane cell parameter of the BaO/SrO superlattice in 
its paraelectric phase is close to the average of the BaO 
and SrO cell parameters (see Supplementary Informa- 
tion). At this relaxed cell parameter, BaO will be under 
compressive strain and SrO under tensile strain. Since 
such epitaxial strains will favor the c-phase in BaO and 
the a-phase in SrO, we can expect a competition between 
the two orientations of the polarization in the superlat- 
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tice. In the following, all the epitaxial strains reported 
are calculated with respect to the relaxed in-plane cell 
parameter of the superlattice in its paraelectric phase. 

In Figure l.e we show the evolution of the A2„ and 
E u modes with respect to epitaxial strain. We find that 
both of them are unstable for a range of epitaxial strains 
from -0.24% to 1.58%; outside this range either A2 U or 
E„ is unstable. This strongly contrasts with the case of 
individual AO compounds, where the A2„ and E„ modes 
are never unstable together, and the regions in which 
they are respectively unstable are separated by a range 
of epitaxial strain without any instability. This feature 
makes the superlattice even more interesting: Whatever 
the epitaxial strain, we expect the ground-state to be 
ferroelectric. 

This is confirmed by the structural relaxation which 
produces the following sequence of ferroelectric ground- 
state for the BaO/SrO superlattice: (i) for r\ <0.07%, 
the ground-state is ferroelectric with polarization along 
the out-of-plane direction (c-phase); (n) for 0.07%< 
i] <1.15%, the ground-state corresponds to a monoclinic 
phase where the polarization has a component in the in- 
plane direction and one in the out-of-plane direction giv- 
ing rise to a small relaxation of the angle between the 
x and z axis (ac-phase); (in) for ij >1.15%, the ground- 
state is ferroelectric with polarization along the in-plane 
direction (a-phase). 

Having demonstrated that ferroelectricity can be strain 
engineered in different alkaline-earth oxides, we now in- 
vestigate the possibility of inducing ferroelectricity in the 
ferromagnetic binary oxide EuO which crystallizes in the 
same rocksalt structure, as a novel route to multiferro- 
ism. Indeed, as for the non-magnetic oxides, we find that 
EuO becomes ferroelectric under epitaxial strain. Impor- 
tantly, its magnetic state remains ferromagnetic through 
the ferroelectric region. The calculated critical epitaxial 
strains for EuO are rf c ~-3.3% and rfe< v ~ 4.2%. These 
critical strains are larger than those of BaO but remain 
experimentally achievable. In Figure l.f we show the cal- 
culated evolution of the spontaneous polarization with re- 
spect to the epitaxial strain. Here, again, the amplitude 
of the polarization can reach sizable values (60 /iC.cm~ 2 
at -5.5%) with out-of-plane or in-plane orientation for 
respectively compressive and tensile epitaxial strain. 

We hope that our theoretical findings will motivate 
experimentalists to confirm our prediction that ferro- 
electricity and related dielectric functionalities can be 
induced in AO compounds of rocksalt structure using 
epitaxial strain. Indeed, it should be possible to tune 
or modify the properties in thin films and superlattices 
by appropriate choice of substrate lattice constant. If 
the epitaxial strain is bigger than the critical value i] c , 
a ferroelectric ground-state with non-zero polarization 
and piezoelectric response should be observed. If the 
experimental strain is substantial, but smaller than the 
critical value, the material will be paraelectric, but its 



TO mode frequency and dielectric constant should be 
strongly strain-dependent as reported in Figure l.d. 

Strain-induced ferroelectricity in simple rocksalt AO 
oxides comparable to that in ABO3 oxides is not only 
an astonishing fundamental finding but also of direct 
interest for practical applications. Although ferroelec- 
tric perovskites are considered as potentially interesting 
high-k gate dielectrics for MOSFET transistors, a signif- 
icant drawback is their small bandgap which results in 
unsuitable band offsets with Si. This problem is directly 
overcome with strained AO compounds. Interestingly, we 
predict that BaO films epitaxially grown on Si (?/=-2.0%) 
should exhibit a ferroelectric ground-state under appro- 
priate electrical boundary conditions [24]. In the context 
of multiferroism, the possibility to induce ferroelectric- 
ity in simple magnetic oxides is particularly attractive; 
previous searches have focussed on more complex com- 
pounds [15, 17]. The case of EuO is particularly inter- 
esting since together with inducing ferroelectricity, the 
epitaxial strain is expected to increase the ferromagnetic 
Curie temperature [25], providing a route to higher tem- 
perature multiferroics. The possibility to further tune 
the properties in superlattices combining different AO 
compounds offers tremendously promising and still un- 
explored possibilities. 

Method 

For the alkaline earth compounds, we performed first- 
principles calculations within the local-density approxi- 
mation to density-functional theory. We used the plane- 
wave based ABINIT [27] software with norm-conserving 
pseudopotentials [26] and a high convergence was ob- 
tained with a plane-wave energy cutoff of 50 Ha and a 
grid of 6x6x6 special k-points to sample the Brillouin 
zone. The vibrational properties, Born effective charges, 
dielectric tensors and piezoelectric constants were cal- 
culated using the density functional perturbation the- 
ory [18]. 

For EuO, we performed the calculation with the VASP 
code [28] and using the GGA PBE functional [29] with 
the DFT+U [30] technique on the / orbitals of Eu (U=7 
eV and J=l eV). Good convergence was obtained with 
an energy cutoff of 600 eV and a 6x6x6 grid of kpoints. 

Work at Liege University was supported by EMMI and 
MaCoMuFi European Project and work at UC Santa 
Barbara was supported by the National Science Foun- 
dation. 
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